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Abstract 

We study limit cycles of nonlinear oscillators described by the equation x + vF{x) + x = with 
F an odd function. Depending on the nonlinearity, this equation may exhibit different number 
of limit cycles. We show that limit cycles correspond to relative extrema of a certain functional. 
Analytical results in the limits v ^ {) and ^ cxo are in agreement with previously known criteria. 
For intermediate v numerical determination of the limit cycles can be obtained. 
PACS numbers: 05.45.-a, 02.30.Hq, 02.30.Xx, 02.60.Lj 
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I. INTRODUCTION 



The study of Lienard's differential equation 



X + h'f{x)x + a; = 



(1) 



or the equivalent form, given earlier by Rayleigh, 



z + uF{z) + z = 



(2) 



goes back to the work of Rayleigh [Q] motivated by his work in acoustics and to E. and H. 
Cartan , van der Pol and Lienard 0] motivated by their work in electrical circuits. This 



type of equation arises directly in numerous applications and others can be reduced to them. 
The problem of studying the number and location of the periodic solutions for polynomial 
F is a particular case of the general two dimensional problem x = P{x,y), y = Q{x,y), for 
polynomial P and Q, which constitutes Hilbert's 16th problem. The problem we address 
here is the determination of the number and position of limit cycles of the equations above 
for even /, or equivalently, in Rayleigh's form for odd F. The two equations have the same 
number of limit cycles, Lienard's form (1) is obtained by taking the derivative of (2) and 
calling z = X. These equations have a unique equilibrium point z = around which their 
periodic solutions will be nested. One of the most studied equations in this class is van der 
Pol's equation which has a unique limit cycle. Lienard established conditions on F which 
guarantee the existence of a single limit cycle. Many results on existence and uniqueness 
have been established, less is known about the number and location of limit cycles which 
do not satisfy Lienard's conditions H, 01- Lins, de Melo and Pugh conjectured that if 
F is a polynomial of degree 2n + 1 or 2n + 2 then there can be at most n limit cycles. This 
conjecture was proved for small z/, that is, for small departures from the hamiltonian 
case. For z/ — the number and position of the limit cycles is given by the real roots of 
a polynomial obtained from Melnikov's function P, ^. In the present case it reduces to 
finding the roots of 




(3) 
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where the integral is performed over a periodic solution Tq of the hamiltonian system z + z = 
0, a circle in phase space. The conjecture has not been proved for arbitrary values of u 
except for special cases. It holds true if F is a polynomial of degree 3 or 4 and if F 



is odd of degree 5 |TT|. These results state the maximum number possible of limit cycles 
but not the exact number. Precise results on the exact number exist for particular cases, 
see for example [O]. In recent work the case u ^ oo has been studied, a criterion to 
determine the number of limit cycles in this asymptotic regime was given and tested in 
several examples |15[. In a different approach, an algorithm to determine the number 



and position of the limit cycles for all values of u and odd F has been formulated. It is non 
perturbative in nature, it is based on finding the roots of a certain sequence of polynomials 
|16| , ITSf. Fewer results [|l3l are known for the generalized equation x + h'f{x)x + g{x) = 0. 

In the present work we study limit cycles of (2), with F{z) and odd function. We show 
that limit cycles correspond to extrema of a certain functional. For small u we recover 



known analytic results namely, (3); for z/ ^ oo we recover the results of [jT4|, [Tsl. For 
intermediate u we must resort to numerical calculations. The approach is a generalization 
of a method developed for other nonlinear problems |]19|, |20|, ^ . This first approach enabled 
us to show that in cases of cubic F with a unique limit cycle, its location can be obtained 
from a minimum principle Here we show that all limit cycles correspond to extrema 

of a certain functional which would allow us, at least numerically, to count the number of 
limit cycles of a given equation. In Section 2 we derive the variational principle, in Section 
3 we obtain analytically criteria for the small and large u regime. In Section 4 we present 
some examples where we see that approximate numerical determination of the limit cycles 
can be obtained far from the asymptotic regimes. Concluding remarks are made in Section 
5. 

II. VARIATIONAL PRINCIPLE 

For odd F, due to the symmetry of equation (2), the limit cycle extends between a 
minimum Zmin = —A and a maximum Zmax = A- Moreover, in phase space, if the point 
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{z, z) belongs to the limit cycle, then the point (— i, —z) also belongs to it. Therefore we 
may consider the positive upper half i; > of the phase plane, where half a period will 
elapse when going from the points (i = 0, z^m) to (i = Q^Zmax)- Then the equation for 
the limit cycle in phase space {z{z),z) can be written as the nonlinear eigenvalue problem, 
VVz + i'F{p) + 2 = 0, with p(±A) = and p > 0. Here we have called p{z) = z{z) and 
the subscript denotes derivative. The eigenvalue is the amplitude A which appears in the 
boundary conditions. It is convenient to define a new variable u = z/A in terms of which 
the equation is 

^p^ + F{p) + Ru = 0,with p(±l)=0, p>0. (4) 

Two parameters appear naturally, R = A/u and S = vA. We may now construct the vari- 
ational principle. Consider the extrema of the functional 0] of two arbitrary functions, 
p{u) with p(±l) = 0, and given by 

/'i {-\P^-F{p))(t){u) du 

R[p{ul 0H] = ^ ^ , 7 (5) 

j_i u(j)[u)du 



/-I \lsp'fu-np)m 



du 

- (6) 



j-i u(j){u) du 

where the second expression is obtained after integration by parts. Variation with respect 
to p{u) at fixed yields the equation 

Here Fp denotes derivative of F with respect to p. Variation with respect to at fixed p 
yields 

^p'^ + F{p) + Ru = Q, (8) 

that is, the equation for the limit cycles. Extrema of R[p, (j)] satisfy both equations (|^) and 
(I). Now notice that equation (0) can be solved for in terms of p. Its solution is 

Finally, replacing in (|^) we obtain the main result, solutions of the equation for the limit 
cycles are extrema of 

/-i</>H {\v{u)FMn)] - F\p{u)]) du 

R[p\ = ext ^ .1 ^ , (10) 

J_i (p[u)udu 



where the extremum is taken over all positive functions p{u) that vanish at the end points 
and (f){u) is the function of p{u) given by @. If we succeed in finding all the extrema, we 
have found all limit cycles. 



III. SMALL AND LARGE i' LIMITS 



A Small 1/ 



Here we shall see that for small u we recover Melnikov's criterion. From the definition of 
the parameters R and S, we see that small or large u, for arbitrary A, corresponds to small 
or large S respectively. For small z/ or equivalently for small S, 

m F,[p{t)] 



/U 



pit) 



dt 



and 



j\C2pF,-F) du 
•St,\^{l-u^)du 



^-ext,i#^^^^^. (11) 



"1 iFp 
-1 2 p 

Let us calculate the first variation 5R of R with respect to p. We obtain 
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S{R[p + 5p\ - R[p\) = — J ^ 6p [pF,, - F,) 



l-RS^-^ 



(12) 



p 

where we called D the integral in the denominator of (pJ]). The term pFpp — Fp does not 
vanish. Then 6R = for arbitrary 6p if 

p{u) =RS{l-u^). 

We know that this is the correct answer from direct integration of the equation. In the small 
u limit the cycle is approximately a solution of the equation ^p^ + Ru = 0, with p(±l) = 
whose solution is what we just obtained, the ellipse p(u) = RS^fl^^v? . This indicates that 



we should use a trial function of the form p{u) = K\/\ — and search for the value of K 
for which R has an extremum. We first notice that we will find some false extrema which 



we must discard. To see this, observe that with this trial function, bp = bK\J\ — v?^ and 
equation (p!2| ) becomes 

SSR =—(!-—) j ^ SKj^ipFpp - Fp)du, 



and the integral will vanish for some K. This is not the desired solution, it is false and 
obtained for not having swept over all possible trial functions. With this in mind, we go 



back to (|TTD. With p{u) = Ayl^^v? as the trial function, equation ([Tl|) is of the form 

hi A) 
S 9{A) 

where we called 

h{A) = j'{lpF,-F)du, 
/•I 1 

9{A) = ^pFpdu. 

It is easy to verify that Ah' — Ag' = —2g, and we obtain 

dR _ 2gh + A{h'g - hg') _ Ah'{A){h - g) 
dA ~ Si^ ~ Si^ ■ 

Extrema of R occur when Ah'{A) = and when h{A) = g{A). The first condition is either 
A = 0, the trivial solution which is always present or when h'{A) = 0. This solution is 
precisely the false solution which we discard. We retain then the solution h{A) = g{A), 
which is simply 

F\p{u)])du = 0, with p{u) = AVl - u^- 
This is exactly condition (Bl) since we have considered only odd F. 



B Large u 



In recent work Lopez et al.[T3|| study limit cycles of Lienard's differential equation (1) 
in the strongly nonlinear regime. Their approach is based in constructing approximate 
solutions to the differential equation. Here we give a heuristic derivation of their result. 
Since 5* appears only in the exponential in the form Exp(S z(u)), we know, from Watson's 



lemma , that when 6*^00 the leading contribution to the integral comes from the points 
where the term in the exponential has an extremum. Here the term in the exponential is 



z{u) 



Pit) 



-dt. 
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Extrema occur where z'{u) = 0, that is where Fp = 0. For large 5* then, R will be given by 

where p is the solution of Fp{p) = and u is the value of u for which p = p on the orbit. 
The extremum is now taken over u. Since R is positive, possible extrema of R will occur at 
M = 1 if F{p) < or at M = — 1 if F{p) > 0. The value of R at the extremum points, if any, 
for large S, are given then by 

R=\F{p,)\, \F{p,)\, etc. 

In order to see which of these correspond to true extrema, hence to limit cycles, we return 
to the differential equation. For large S, limit cycles are approximately given by, 

F{p) + Ru = (13) 

the solution of which has to be matched to a thin boundary layer. Now that we know the 
possible values of R, we can read directly which are the limit cycles. This is best seen 
in a plot in phase space of ([T3|). Suppose there are several points p at which Fp = 0. 
Since F is odd, take the positive values and label them in growing order, pi < P2 < Ps---- 
Correspondingly we have different values of R. Since limit cycles must be nested, and, for 
the systems considered, the derivative dp/du can vanish only at one point, the only limit 
cycles will be those for which > \F{pk)\ for all k < j. This is best seen in Fig. 1 

where we show a case with two allowed limit cycles. In the asymptotic regime when u ±1 
the horizontal coordinate tends to R = F{p). The trajectory of each limit cycle is indicated 
by the arrows. From the definition of R we know then that in the limit z/ ^ oo, for each 
allowed cycle, the amplitude grows as Ai = \F{pi)\i'. Moreover we can read the maximum 
value of p in each cycle. The maximum p are solutions of F{p) + = 0. The values of 

p determined in this way corresponds to the amplitude of the associated Lienard equation 
(1). Thus we have recovered the solution of [|1^, |15], which gives support to the conjecture 
of Lins, de Melo and Pugh. 
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IV. EXAMPLES 



Having seen that we recover the small v and large v limits, we now give numerical results 
for arbitrary values of v in simple examples. Here too, false extrema may appear due to the 
impossibility of sweeping over all trial functions. This becomes evident, as in the z/ — > 
limit, by considering the first variation of R with respect to p, we find that for arbitrary S", 
R[j) + — -R(p) vanishes when 



1 

du 5p 



F F 
p p2 



Rt-{-p{t)FMt))-W)) 



0. 



-1 

As in the small u limit, the term in square brackets does not vanish, so for arbitrary 6p the 
first variation vanishes when the term in curly brackets vanishes. Taking the derivative of this 
term one can see that it corresponds to the equation for the limit cycles. In practice though, 
when sweeping over a restricted variety of trial functions the integral may vanish at other 
points, which calls for some care. To illustrate the use of the variational principle we shall use 



for all u a simple trial function with only one parameter, the ellipse p{u) = K\J\ — u^. For 
each value of 5* we insert this trial function in (p!OD and sweep in K to find all the extrema. 
We begin with small S where we identify the true minima and follow their evolution as 5* is 
increased. We obtain a table of extrema i?2, -R3, •••) for each S", from which we compute 



the values z/j = ^SjRi and the corresponding amplitude Ai = y/SRi. As we show below 
with this very simple trial function, appropriate for z/ — > 0, we obtain the approximate 
location of the limit cycles for u well beyond the region of validity of this trial function. 
As a first example we take 

r./ N 4 43 85 
F{p) = ^P- 3P + • 

For 1/ ^ 0, predicts the existence of two limit cycles, of amplitudes A = 1 and A = 2. 
As V increases the amplitudes change. In Fig. 2 we show the amplitudes as functions of v. 
The continuous line is the amplitude obtained from the direct integration of the differential 
equation (Q). Unstable limit cycles are obtained by integrating the equation with u < Q 
since this only changes the stability properties of the cycles. The origin is a stable fixed 
point, the inner limit cycle unstable and the outer stable. The dots indicate the values of 
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the amplitude obtained following the extrema of R[p]. The agreement is close, in spite of 
having used a one parameter trial function. 
A more interesting example is provided by 

Here (3) again predicts the existence of two limit cycles for z/ of amplitudes A = 



41 - l)/5 ^ 1.039 and A = y(V4T+ l)/5 ^ 1.217. For this equation it is known that 
for large u there is no limit cycle. A bound on the value of u for which no limit cycles 



exist is known ||2^. The criterion for z/ ^ oo also indicates that no limit cycles exist in this 
regime. Numerical integration of the differential equation shows that, as u increases these 
two limit cycles merge and disappear. In Fig. 3 we show the amplitude as a function of u. 
The continuous line is obtained by direct numerical integration and the dots were obtained 
variationally. There is qualitative agreement, the variational result predicts the merging of 
the limit cycles. The exact cycles are very different from the circles we use as trial function 
so this trial function is not adequate at larger i/. 



V. SUMMARY 

We have shown that all limit cycles of (0) correspond to extrema of a certain functional. 
The exact position of each limit cycle is obtained when the trial function coincides with the 
solution, otherwise an approximate estimation can be made. From the variational expression 
analytical results can be obtained in the two asymptotic limits, z/ — and u oo. In these 
asymptotic regimes we obtain both the amplitudes of the Rayleigh and Lienard form of 
the equations. In the small u limit we reobtain the known criterion, namely Melnikov's 
integral. In the large u regime our results coincide with that obtained in recent work. In the 
intermediate regime the number and position of limit cycles can be obtained numerically. 
Even with simple trial functions relatively close estimates are obtained. Whether it is 
possible to count analytically the number of extrema, hence of limit cycles, is a question 
that remains to be answered. 
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FIG. 1: Graph of allowed limit cycles in the u ^ oo regime. 
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FIG. 2: Amplitude ofthe limit cycles for = (4/5)p- (4/3)p=^+(8/25)p^ The continuous 
line is the result of numerical integration of the differential equation. The dots were obtained 
variationally with a simple trial function. 
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FIG. 3: As in Fig.2, for F{p) =p - ^41/9p3 + pS. 
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